Effects of mesonic correlations in the QCD phase transition 
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O ' The finite temperature phase transition of strongly interacting matter is studied 

o : 

within a nonlocal chiral quark model of the NJL type coupled to a Polyakov loop. In 
contrast to previous investigations which were restricted to the mean-field approx- 
imation, mesonic correlations are included by evaluating the quark-antiquark ring 
sum. For physical pion masses, we find that the pions dominate the pressure below 
the phase transition, whereas above T c the pressure is well described by the mean- 

JL , field approximation result. For large pion masses, as realized in lattice simulations, 

M-t, 

r- ' the meson effects are suppressed. 

rn 

> ■ The interest in the phase diagram of quantum chrome-dynamics (QCD) has received new 

: 

impulses from the recent results of heavy-ion collision experiments at RHIC Brookhaven [1], 

m 

revealing that hot hadronic matter at temperatures not too far above the critical tempera- 
ture T c behaves like a perfect liquid, rather than a weakly coupled plasma. A theoretical 
description of these investigations therefore requires a non-perturbative approach, which also 



X 



provides a proper understanding of the chiral quark dynamics and the confinement mecha- 
nism. Until now, the only method which is directly based on QCD and which meets these 
requirements is lattice gauge theory. Unfortunately, the application of lattice results to ex- 
perimental data is complicated by the fact that most lattice calculations are performed with 
rather large quark masses, leading to unphysically large pion masses. Therefore, in spite 
of recent progress in this aspect, additional methods are required to extrapolate the lattice 
results to the physical masses. To bridge the gap there are major efforts within chiral per- 
turbation theory (xPT) (see, e.g., [2]). However, %PT is obviously not suited to describe 
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the region of the phase boundary or even beyond. In this regime, one has to employ other 
approaches for the chiral extrapolation, like QCD Dyson-Schwinger equations (see, e.g., [3]) 
or effective models which share the relevant degrees of freedom with nonperturbative QCD 
[4]. These approaches have the additional advantage that they can straightforwardly be 
extended to nonzero chemical potential, which is still a serious problem on the lattice. 

In this article we want to discuss an effective model of low-energy QCD, capable of 
describing the chiral as well as the deconfinement transitions. As a basis we use the PNJL 
model [5-8] which generalizes the well-known Nambu-Jona-Lasinio (NJL) model [9] for the 
chiral quark dynamics by coupling it to the Polyakov loop, which serves as an order parameter 
of the deconfinement transition. To a large extent, this removes one of the most disturbing 
features of the original NJL model, namely the pressure contribution of unconfined quarks 
in the hadronic phase. In spite of the simplicity of the model a remarkable agreement with 
the results of lattice QCD thermodynamics [10, 11] has been obtained [5]. 

However, this comparison was not entirely consistent: Whereas unphysically large values 
of the current quark masses have been used in the lattice simulations [10, 11], physical va- 
lues have been employed in the PNJL model analysis of Ref. [5]. Moreover, after successfully 
removing (most of) the unphysical quark degrees of freedom from the confined phase, the 
PNJL model treated in mean-field approximation as in Ref. [5] does not contain any degree 
of freedom in this regime. Obviously, this is a rather poor description of the hadronic phase 
at finite temperature where mesons are expected to become relevant. Hence, the good agree- 
ment of the PNJL results with the lattice data could partially be accidental in that mesonic 
correlations have been neglected in the PNJL analysis, while in the lattice calculations they 
are suppressed by the large current quark masses. To get a consistent picture it is thus 
important to go beyond the mean-field approximation and to include mesonic correlations. 
This can be done systematically within the framework of a 1/N C expansion scheme [12], 
which has successfully been applied to the NJL model, e.g., in Refs. [13]. 

In the present work we suggest a l/N c improvement of the PNJL model which is necessary 
to disentangle hadronic contributions in the vicinity of and below the chiral/deconfinement 
phase transition. Moreover, the existence of bound states above T c may be crucial for the 
understanding of the properties of the strongly coupled quark-gluon plasma [14]. A consis- 
tent inclusion of hadron gas contributions should thus include the dissociation of hadrons 
being bound states of quarks and antiquarks below the transition into resonant continuum 
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correlations somewhere above it. This hadronic Mott-transition [15] has been discussed first 
in Refs. [16] and was formulated within the NJL model in Ref. [17] employing the concept of 
spectral functions and in-medium scattering phase shifts within a Beth-Uhlenbeck approach, 
see [18]. 

In this letter, we restrict ourselves to the case of zero chemical potential. The central 
quantity for our analysis is then the thermodynamic potential at finite temperature which can 
be written as a sum of a mean-field part and a part which describes the mesonic correlations, 

n(T) = n mi (T) + n COTT (T)-n . (1) 

As usual, we have introduced an irrelevant constant f2 , which is chosen such that f2(0) = 0. 
Technically, the correlation part f2 corr (T) corresponds to a ring-sum of quark-antiquark loops 
(see [17, 19]). For the evaluation of (1), we consider a nonlocal generalization of the PNJL 
model. The nonlocal four-point interaction is chosen in a separable form motivated by the 
instanton liquid model (ILM) [20] where it results from the internal nonlocal structure of 
the nonperturbative QCD vacuum. In the ILM, the nonlocality is represented by the profile 
function of the quark zero-mode in the instanton field and depends on the gauge. We use 
a Gaussian ansatz as one of the simplest functional forms of the nonlocality which has a 
similar behavior as the zero-mode profile obtained in a gauge invariant manner [21]. This 
choice guarantees convergence at all orders without an additional regularization procedure. 
The quark sector of the nonlocal chiral quark model is described by the Lagrangian 

C q = q{x){iip - m c )q{x) + j[J 2 a (x) + J 2 {x)\ , (2) 

where m c is the current quark mass, and = <9 M — iA^ the covariant derivative with a 
background gluon field = A ^ = S^ A . The nonlocal quark currents are 

Ji{x) =J d 4 x 1 d 4 x 2 f(x 1 )f(x 2 ) q(x - x ± ) r 7 q(x + x 2 ), (3) 

with the vertices Y a = 1 and = ry 5 r a , a = 1, 2, 3. 

After linearization of the four-fermion vertices by introducing auxiliary scalar (a) and 
pseudoscalar (n a ) meson fields the quark sector is described by the Lagrangian 

2 ~ 2 

A- = - m c )q{x) - + J*(x)a(x) + n a {x) J$(x) . (4) 

To proceed, we single out the nonzero mean-field value of the scalar field by the decomposition 
a = (7 + er m f so that n a and a denote only the fluctuating parts of the fields ((7r a ) = (a) = 0) 
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describing mesonic correlations. The scalar mean field gives a dynamical contribution to the 
quark mass which in Euclidean momentum space is M{p 2 ) = m c + rridf 2 (p 2 ) with / 2 (p 2 ) = 
exp(— p 2 / A 2 ) being the (Fourier transformed) Gaussian form factor in Euclidean space. The 
amplitude = — <7 m f is an order parameter for dynamical chiral symmetry breaking. The 
model parameters of the quark sector are fixed by using the pion mass and the weak pion 
decay constant and choosing a chiral condensate value in accordance with the limits from 
QCD sum rule analyses. The chiral condensate is obtained from the non-perturbative part 
of the quark propagator, S np (p) — + M(p 2 )) -1 — {j> + ^ c ) _1 , i.e., after subtracting the 
perturbative part. For details see, e.g., [22, 23]. For the analyses in the present paper, we 
employ a parameter set from Ref. [22] for — (qq) 1 ^ 3 = 240 MeV given by m c = 5.8 MeV, 
A = 902.4 MeV and GA 2 = 15.82. 

The mean-field thermodynamic potential reads 

anfCO = -4 E / iog[K) 2 + ^ 2 + M 2 (K) 2 )] +!3i+w(<M) , (5) 

i=0,± , 
' k,n 

where the notation J kn = TJ2 n f d 3 k/(2n) 3 has been introduced. <3> denotes the Polyakov 
loop variable, which is given by $ = |Tr c e l< ^ T . Here = A A = iA is related to the 
(Euclidean) background gauge field. In Polyakov gauge it is diagonal in color space, i.e., 
= 3 A 3 + 08-^8- Following Ref. [8], we require $ = $ to be real with real 03,08- As a 
consequence 08 = and we are left with one variable 03 [8]. 

Due to the coupling to the Polyakov loop the fermionic Matsubara frequencies u n = 
(2n+ l)vrT are shifted: 



CO 



± 



LO n ± 03, CJ° = 0J n . (6) 



For the Polyakov loop potential U(&, $) we adopt the the logarithmic form of Ref. [8], which 
has been fitted to the quenched lattice data of Ref. [24]. 

The order parameters (mean field values for and 03) are obtained by minimization of 
the mean-field part of the thermodynamic potential = 0, = 0. 

In Fig. 1 we show the resulting temperature dependence of the quark condensate {qq) T 
(normalized to its vacuum value) and of the Polyakov loop expectation value. We also show 
the corresponding results for the case of uncoupled (pure) quark and gluon sectors. In the 
pure quark model the critical temperature for chiral restoration is T c = 116 MeV, whereas 
the pure gauge sector has a critical temperature for deconfinement Td = 270 MeV, fixed 
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from lattice data for the Polyakov loop. When the quark and gluon sectors are coupled 
these temperatures get synchronized so that T c « « 200 MeV. 




T (MeV) 

Figure 1. Quark condensate (normalized to its vacuum value) in the nonlocal NJL model (green 
dotted line) and the nonlocal PNJL model (red solid line) compared to the Polyakov loop in pure 
gauge theory (blue dash-dotted line) and in the nonlocal PNJL model (black dashed line). 

We now include fl C orr to study the effect of mesonic correlations. The properties of mesons 
in the PNJL model have been studied in Ref. [7] and the generalization to our nonlocal model 
is straightforward. Key ingredients are the quark-antiquark polarization loops Hm- Using 
the notations D(k 2 ) = k 2 + M 2 (k 2 ), k % n+ = (u l n + u m , k + p), and k\ = (cj* , k) they can be 
written as 

rV(p> m ) = m fEj kn ^^^^j K + -K±M((K + ) 2 )M((ki l ) 2 ))]. (7) 

The mesonic contributions to the thermodynamic potential are then (to lowest order) given 
by the ring sum [17], 

fWT) = IT / Ml-Gn M (p> ro )] (8) 

with the mesonic degeneracy factor du- 

The fermionic Matsubara sum, the three-momentum integration and the sum over the 
three color modes i = 0, +, — in (7) are evaluated numerically with a finite result for 11^ (p, z/ m ) 
to be inserted in Eq. (8). The bosonic Matsubara sum in (8) and the three momentum 
integral are also performed numerically. 

Our model predictions for the the pressure P(T) = — fi(T), divided by the Stefan- 
Boltzmann limit, are displayed in Fig. 2. For comparison with the full result we also show 
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the mean-field result and the mean-field plus pion contribution as well as the result for an 
ideal pion and sigma gas with the masses fixed at their vacuum values. We find that at low 
temperatures the mean-field (i.e., quark) contribution is suppressed and the pressure can 
be well described by a free pion gas. Near the critical temperature the a meson gives an 
additional visible contribution whereas already after T > 1.5 T c the mesonic contributions 
are negligible and the quark-gluon mean-field dominates the pressure. 
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Figure 2. Scaled pressure p/psB in the nonlocal PNJL model with the physical pion mass: mean 
field contribution (green dotted line), mean field + pion (blue dashed line), mean field + pion + 
sigma (black solid line). The red dash-dotted line denotes the scaled pressure of an ideal pion + 

sigma gas with fixed masses. 
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Figure 3. Scaled pressure p/psb in the nonlocal PNJL model with the physical pion mass (thick 
lines) and with m c /T — 0.4 (thin lines): mean field contribution (green dotted line), mean field + 

pion + sigma (black solid line). 

So far we have used the parameter set of Ref. [22] with m c = 5.8 MeV, corresponding to 
the physical pion mass of 140 MeV in vacuum. On the other hand, in most lattice calculations 
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the masses are considerably larger. Thus in order to perform a meaningful comparison with 
lattice results, we should repeat our model calculations using quark masses similar to the 
lattice ones. To be specific, we choose the current quark mass to scale with the temperature 
as m c = 0.4 T, mimicking the situation in the lattice calculation of Ref. [25], where the up 
and down quark masses behave in the same way. 

Our results for the scaled pressure are displayed in Fig. 3. For comparison we show 
again the results obtained with m c = 5.8 MeV (thick lines). With m c = 0.4 T (thin lines) 
the qualitative behavior remains unchanged: Above T c the total pressure (black solid line) 
converges quickly to the mean-field result (green dotted line) , whereas below T c the pressure 
is dominated by the mesonic contribution. Quantitatively, the meson contributions are of 
course strongly suppressed in the case of the heavier quark masses. Also note that the mean- 
field result is affected as well. In particular T c rises from 202 MeV for m c = 5.8 MeV to 
231 MeV for m c = 0.4 T 
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Figure 4. Scaled pressure p/psB as a function of T/T c : non-local PNJL model with physical pion 
mass (black solid line) and with m c /T — 0.4 (blue dotted line). Red dashed line: Lattice data for 

two-flavor QCD with staggered quarks [25]. The shaded region is an estimated continuum 
extrapolation of these data for massless QCD [25]. Points: Lattice data for two-flavor QCD with 
Wilson- type quarks [11] for N t = 6 (open symbols) and JV* = 4 (green filled symbols). The data for 
the pressure [11] have been divided by the Stefan-Boltzmann limit for N t = 6 and Nt = 4, 

respectively, as given in Ref. [11]. 

In Fig. 4 we compare our model results for the scaled pressure with lattice simulations. 
Our model results are indicated by the black solid line and the blue dotted line, which 
correspond to the physical parameterization (m c = 5.8 MeV) and the case m c = 0.4 T, 
respectively. Both curves include mesonic correlations and are identical to the black solid 
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lines in Fig. 3. Note, however, that they are now displayed as functions of T/T c , with T c 
being different for the two parameterizations. 

The red dashed line indicates the result of Ref. [25] obtained on a 16 3 x 4 lattice with 
improved staggered fermion actions for two flavors. As mentioned above, in this calculation 
the quark mass is unphysically large and scales with the temperature, m Ui d/T = 0.4. At 
T > 2T C our results are systematically above the lattice data. Since the latter are obtained 
with four time slices only (N t = 4), this is probably a lattice artifact. This interpretation 
is supported by the lattice results of Ref. [11] for Wilson fermions with N t — 4 (green filled 
symbols) and N t — 6 (open symbols). While the former almost coincide with the results of 
Ref. [25], the latter indicate higher pressures at large T. Although the Wilson data have 
been extracted on "lines of constant physics", i.e., for T-independent quark masses and are 
therefore not entirely comparable to the staggered fermion data, this tendency should be 
correct. Indeed, the authors of Ref. [25] tried a continuum extrapolation for massless QCD 
(shaded area), which is in fair agreement with our model results for both, physical and large 
quark masses. 

Whereas in our model, even for m c = 0.4 T, there remains a visible pion contribution 
to the pressure down to T/T c m 0.2, the lattice pressure of [25] vanishes at T/T c = 0.6. 
This discrepancy has a trivial explanation by the fact that the lattice pressure has been 
obtained by the "integral method" which leaves one integration constant undetermined. In 
Ref. [25], this constant has been fixed by the choice that the pressure vanishes at T/T c = 0.6. 
On the other hand, chiral perturbation theory predicts that at very low temperature the 
pressure is well described by an ideal pion gas, in good agreement with our model. Hence, 
if the integration constant on the lattice had been fitted to xPT ■, rather than setting the 
pressure at T = 0.6 T c equal to zero, they would be in good agreement with our results at 
this point. On the other hand, our approach underestimates the lattice data in the region 
0.9 < T/T c < 1.6. This may be attributed to our neglect of hadronic resonances, other than 
pion and sigma. 

In summary, we have extended previous studies of the thermodynamics of low-energy 
QCD within Polyakov-loop NJL models to include mesonic correlations. To that end, we 
have evaluated the pion and sigma contributions to the pressure by calculating the ring 
sum. We have thereby employed a nonlocal four-fermion interaction of the instanton liquid 
type. The pionic contribution dominates the pressure in the low temperature region. In 
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this regime, the pressure is quite insensitive to the details of the interaction and agrees 
almost exactly with that of an ideal pion gas. At temperatures > T c , on the other hand, 
the mesonic contributions die out. For unphysical pion masses, as obtained in most present 
lattice calculations, mesonic correlations play only a minor role, even below T c , in accordance 
with the lattice results. 

For the future we plan to include the back reaction effect of mesonic correlations to the 
equation of motion of mean fields. This can lead to a lowering of T c . We also plan to study 
nonzero chemical potentials, which may require a nontrivial extrapolation of the Polyakov- 
loop potential into this regime [26]. We then should include baryonic degrees of freedom as 
well. 
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3<D<DEKT ME30HHHX KOPEJIJI51Ii;HH B OA30B0M 

nEPEXO^E KXA 

JJ,. BjianiKe, M. By6ajia, A.E. Pa,a;:>Ka6oB, M.K. Bojikob 

<J>a30BbiH nepexo,!]; CHjibH0B3aHM0fleHCTByK>in,eii MaTeproi npii kohchhoh TeMnepaType 
H3yHeH b HejiOKajiBHoft KHpajibHoii KBapKOBoii Mo^ejiH Tnna HHJ1 coeflHHeHHoii c neT- 
jie-H IlojiHKOBa. B OTjiiiniie ot cymecTByiomnx nccjieflOBamiH b paMKax cpe^Hero nojia 

B MOflejIB BKJIIOTeHbl Me30HHbie KOppejIHUHH C nOMOIUblO KBapK-aHTHKBapKOBOH KOJIb- 

u;eBoft cyMMM. Jlflft (pH3HHecKiix Mace nnoHa HaMH o6Hapy>KeHO, hto ihiohbi ,a;aiOT ,n,o- 
MHHHpyromHH BKjia/i; b ,a;aBjieHHe, Tor^a KaK Bbiine T c flaBjieHiie xopoino oniicbiBaercH 
pe3yjibTaTOM BbiHiicjiemiH b npn6jiHJKeHHH cpe,a;Hero nojia. J^jisi 6ojibinofi Maccbi niiOHa, 
KaK b BbiHHCJieHHHx Ha peineTKe, Me30HHbie scpcpeKTbi no^aBjieHbi. 



